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Abstract 



The chaotic spike train of a homochnic dynamical system is self-synchronized 
I by re-inserting a small fraction of the delayed output. Due to the sensitive 

o ■ 

' nature of the homoclinic chaos to external perturbations, stabilization of very 

o ■ 

long periodic orbits is possible. On these orbits, the dynamics appears chaotic 
over a finite time, but then it repeats with a recurrence time that is slightly 
^ ' longer than the delay time. The effect, called delayed self-synchronization 

' (DSS), displays analogies with neurodynamic events which occur in the build- 

' up of long term memories 



PACS numbers : 05.45.Xt, 05.45.Vx, 42.65.Sf. 
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The last decade has seen a great interest in controlhng chaos using small perturbations 
||1|]. It started with the seminal paper by Grebogi, Ott, and Yorke in which they proposed a 
method of stabilizing unstable periodic orbits using tiny, yet carefully chosen perturbations 
to an accessible system parameter. Subsequently, other methods were proposed, as the 
delayed feedback due to Pyragas P|, who showed that re-inserting a time-delayed version 
of the output back into the system can at certain conditions stabilize some of its periodic 
orbits. However, that method does not provide a reliable procedure that guarantees the 
stability of a chosen orbit, in particular it is exceedingly difficult to stabilize long periodic 
orbits. An improved version of the delayed feedback was represented by the adaptive control 

experimentally implemented via a Taylor expansion which consists in feeding a delayed 
fraction of the output as well as its variation rate ||^. 

In this paper we show that for a chaotic system displaying a continuous return to a saddle 
focus (Shilnikov chaos [0) and a high sensitivity to external perturbations applied near 
that focus, a long-delayed feedback can indeed be used to stabilize very complex sequences 
of pulses. We demonstrate this in numerical simulations and experiments with a C02 
laser as well as on a return map model of a chaotic pulse generator. However the validity 
of this scenario is much broader, including possible neurodynamic implications. We call 
this effect "delayed self-synchronization" (DSS), insofar as the time signal appears as a 
train of N erratically distributed spikes (N being the ratio between the delay time and the 
average inter-spike interval) which repeat themselves with the same interspike intervals. The 
DSS phenomenon should by no means be confused with the synchronization of two distinct 
systems 0. 

Standard control methods are of limited effectiveness for long times insofar as their 
complexity increases with the length of the periodic orbit, as shown e.g. by the Hunt's 
occasional proportional feedback applied to long period orbits g]. On the contrary, in DSS 
each return to the saddle focus permits an independent control of the corresponding inter- 
spike interval, and such an operation can be updated as long as one wishes, without any 
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increase in complexity. In our laboratory implementation (Fig.l), a single- mode CO2 laser 
with intensity feedback is tuned to the parameter range yielding homoclinic chaos ||^. The 
intensity output of the laser consists of a sequence of almost identical spikes, repeating at 
erratic times because of the homoclinic scrambling, with an average period of 500 fis. (Fig. 
2a). Through a delay unit described elsewhere a small delayed perturbation (a few 
percent of the intensity output) is added to the feedback signal responsible for homoclinic 
chaos. As figures 2-b,c show, this small feedback creates a sequence of spikes, which is 
periodic with a period slightly larger than the delay time Td of the feedback loop (DSS). 
In fact, the DSS adjusts in such a way that the spikes of the delayed signal arrive through 
the delay line in the time interval of largest susceptibility, which occurs around r ~ 150/is. 
before the next large spike whence = + t . The extra time or "refractory time" r, 
which in fact corresponds to the duration of the quenched (zero intensity) time interval for 
each spike, is also measured by the width of the correlation function of the free running laser 
(Fig.3-a). 

The autocorrelation function of the DSS signal has large revival peaks separated by the 
peiod Tj. (Fig. 3-b). Figure 4 shows that the minimal signal e necessary for DSS is constant 
for long delays, down to a delay of r ~ 150/is. For delays below the intrinsic refractory time 
r, the DSS threshold increases dramatically. 

DSS can be further illustrated by the space-time representation (STR) introduced else- 



where |]TT|. In this representation, the time series is split into pieces of length T^, which 
then are stacked together as different "snapshots" of a ID spatiotemporal system. Thus, 
every line of this representation is mapped onto the next line. The STR for the free-running 
system and for the DSS regime are shown in Fig. 5. When DSS is achieved, the STR of the 
dynamics changes drastically. 

We also investigate DSS by a numerical experiment with the model of the homoclinic 



chaos in CO2 laser, introduced earlier [12 
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Xi = koXi{x2 — 1 — ki sin^(x6)) 

X2 = -11X2 - 2kQX2Xi + gx'i + xa + pq 

X3 = -71X3 + X5 + gx2 + po 

xa = -12XA + gx^ + -23:2 + zpQ 

^5 = -722:5 + ZX3 + gxA + zpo 

xg = -/5(x6 + &0 - rXi) 

Here Xi is the laser intensity, X2 is the gain, X3 , X4, X5 are variables of the gain medium, 
x6 is the feedback voltage, and Pq is the pump rate. The delayed feedback is introduced in 
the equation for the feedback voltage via Xi = xi{t) + e ■ {xi {t — Td) — xi), where xi is 
the mean value of xi{t). Parameters of the model at which the regime of homoclinic chaos 
is observed are k^ = 28.57, k^ = 4.55; 71 = 10.06, 72 = 1.06, g = 0.05, po = 0.016, z = 10, 
r = 70, bo = 0.173, (3 = 0.42. 

In numerical simulations, we are not limited in the length of the time delay and the time 
series. We ran very long simulations (up to 5-10^ non-dimensional time units) with long time 
delays up to 8000. The numerical time series (Fig. 6) are very similar to the experimental 
ones. The auto-correlation function also exhibits sharp peaks separated by a time interval 
Tr that is slightly larger than T^. The magnitude of the peaks very slowly decays with time, 
which indicates the slow evolution of the periodic orbit. 

Notice that, once the dynamical system has become periodic with period Tr = + r, 
that is, xi(t) = xi{t — Tr), than by Taylor expansion we can write Xi = {1 + e) ■ xi{t) + e ■ 
(^)r — e ■ Xi. Thus for £ = — 1, the delay feedback effect disappears as zero-th order term, 
as expected from Pyragas , but a first order correction proportional to the time derivative 
of xi represents the adaptive correction of Ref. 0] in the approximation of Ref. 0. 

The DSS mechanism can be elucidated using the following simplified model. Without 
time-delayed feedback, after a pulse is produced, the phase trajectory spends some amount 
of time (which fluctuates because of the homoclinic chaos) near the homoclinic saddle before 
it leaves its neighborhood and produces the next pulse. A close inspection of the time series. 
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both experimental and numerical (see Figs. 2-b,c and 6), indicates that the arrival of a 
pulse through the delay loop triggers the escape of the phase trajectory of the laser system 
from the neighborhood of the homoclinic saddle. Thus, the following simplified model can 
be proposed. Consider a system that generates a pulse at time ti , which is determined by 
a map 

ti = ti_i + / (t,_i - t,„2) (2) 

If the function / (•) is chosen such that the map Xn+i = f produces chaos, then the 
sequence of pulses will be chaotic. Such system has actually been implemented in electronic 
circuitry Let us assume now that the pulses are inserted back into the system after 



a certain time delay Tq, so if a delayed pulse arrives at a time tp between +5 and t,, 
then it triggers a new pulse at tp instead of ti, and the system is reset so that new ti = tp. 
We introduce the small offset S in order to avoid spurious generation of pulses with very 
small time distance, which would effectively destroy pseudo-chaotic pulsation. In fact in the 
experiment this "refractory time" was — = 150//S. We iterated this system numerically 
using the logistic map / (x) = ax {1 — x) with a = 3.8 and found that the system typically 
settles on a periodic orbit with the period Tq as shown by its autocorrelation function (Fig. 7). 
The particular orbit will obviously depend on initial conditions, in particular, orbits with 
periods much smaller than Tq, can also be stabilized, however they are not typical. Note, 
in this case the strength of the feedback is irrelevant as long as the magnitude of the pulses 
passed through the delay line is sufficient to trigger the next pulse. This agrees with the 
experimental observation that the threshold of synchronization is independent of Tq for large 
enough Tq (Fig. 4). We also checked the influence of noise on DSS. If noise was added to 
the sequence of the pulses at the levels below the threshold for pulse generation, it would 
not affect DSS. However, if a small white Gaussian noise is added to the inter-pulse interval, 
the revival peaks eventually decay to zero, albeit very slowly. 

We have thus introduced a conceptual model for the phenomenon of delayed self- 
synchronization which leads to stabilization of very long-periodic orbits in a chaotic system. 



Such long orbits are known as a signature of the so-called pseudo-chaos observed in systems 



with discretized state space [0 . In those systems, the short-term behavior is indistinguish- 
able from a real chaos, however the orbit returns to the initial point after some time and 
then repeats itself. Obviously, in the systems considered here, the state space is continuous, 
so the analogy between pseudo-chaos and DSS cannot be carried too far. 

The ability to synchronize very long and complex periodic orbits is of particular interest 
in relation to the recent discovery of the neuronal mechanism of transformation of short-term 
memories into permanent (long-term) memories via so called synaptic reentry reinforcement 
(SRR) [|15|. In those experiments, it has been demonstrated that neurons responsible for 
learning, repeatedly "replay" certain spiking patterns in order to establish permanent con- 
nections (synapses) among neurons. We believe that stabilization of the complex periodic 
orbits in simple chaotic systems presented in this paper can serve as a paradigm for SRR in 
more complex neural systems. In our forthcoming work we will consider stabilization of long 
periodic orbits in the Hindmarsh-Rose neural model We believe that the homoclinic 
chaos observed in that system as well as in real biological neurons |T^, is analogous to the 
laser chaos described here, and thus is should be susceptible to the phenomenon of delayed 
self-synchronization. 

Authors are indebted to S. Boccaletti and N.F.Rulkov for useful discussions. F.T.A., 
R.M. and E.A. acknowledge partial support from the European Contract No. HPRN-CT- 
2000-158. A.D.G. is supported by European Contract No. PSS 1043. L.T. wants to thank 
Istituto Nazionale di Ottica Applicata (Florence) for warm hospitality. The work of L.T. 
was partially supported by the Engineering Research program of the Office of Basic Energy 
Sciences at the U.S. Department of Energy, grants DE-FG03-95ER14516 and DE-FG03- 
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Figures 




Figure 1 

FIG. 1. Experimental setup consisting of a CO2 laser with a feedback loop (dashed line) impos- 
ing a regime of homoclinic chaos, and a delay unit. EOM, electro-optic modulator; Det, HgCdTe 
detector; r and hQ, gain and bias of the amplifier in the feedback loop; e and Delay, coupling factor 
and delay units. 
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figure 2 

FIG. 2. A sequence of homoclinic spikes in the output intensity of a CO2 laser in the free 
running regime (a) and DSS for two different delays (1 and 4 ms, (b) and (c) respectively). A 
thick arrow denote the zero intensity level. The two horizontal bars for 1 and 4 ms show the role 
of the refractory time r = IbOfis. The shorter bar is the imposed delay (T^), the larger one, with 
r added, is the effective delay {Tr) that characterizes DSS. 
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figure 3 



6 8 10 12 14 16 18 
time (ms) 



FIG. 3. Autocorrelation function of the laser intensity for the free running case (a) and for DSS 
with Td = 4ms (b) . The recurrence time for the revival of the correlation is the sum of plus 
the refractory time r. 
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figure 4 ^'^^^ 
FIG. 4. Minimal DSS signal (in percent of the peak to peak output) versus the delay time. 
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figure 5 

FIG. 5. Space-time plots of the transition from free running to DSS regime with a delay time 
T^i = 4ms (indicated by the arrow). The space cell corresponds to a single recurrence time T^, the 
time coordinate is a sequence of integer values corresponding to the successive delay units. 
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figure 6 

FIG. 6. Numerical simulations of the laser equations (1) with time delay = 2ms: solid line, 
a section of the time series of the laser output; dotted line - the same time series shifted forward by 
the delay time and displaced vertically by —1. As clearly seen, the spikes of the delayed signal 
give rise to the rapid escape of the laser intensity from the floor at 0.5 (denoted by a thich arrow), 
corresponding to the saddle-focus. 
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figure 7 

FIG. 7. Autocorrelation function of the chaotic pulse generator (2) with delayed feedback. 
Sharp spikes are separated by the delay time Tq = 1000. 
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